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Abstract 

The properties of the D*j{2317)~^ and D^5j(2460)"'" mesons are studied using 125 fb^^ of e"'"e~ cc 
data cohected by the BaBar experiment. Preliminary mass estimates of [2318.9 it 0.3 (stat.) it 
0.9 (syst.)] MeV/c^ and [2459.4 ± 0.3 (stat.) ± 1.0 (syst.)] MeV/c^ are obtained. Searches are 
performed for the decay to the Df meson along with one or more vr*^, vr''", or 7 particles. A search 
is also performed for neutral and doubly-charged partners. 
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1 INTRODUCTION 



The -D*j(2317)"^ meson, discovered by this collaboration [1] and confirmed by others [2, 3, 4], 
and the Z)sj(2460)''' meson, observed by the CLEO collaboration [2] and confirmed by this col- 
laboration [5] and others [3], has reawakened interest in the study of the spectroscopy of charm 
mesons. 

Presented in this paper is an updated, preliminary analysis of these states using 125 fh~^ of 
e~^e~ cc data collected by the BABAR experiment. From this analysis we present new estimates 
of the D*j{2317)'^ and Dsj(2460)^ masses and the branching ratios of L>sj(2460)"'" decays to Dfj 
and Z)+7r+7r~ with respect to its decay to Dfir'^j. In addition, we search for new decays of either 
meson involving a meson accompanied by 7r° and tt^ mesons and photons. An estimate of the 
-Dsi(2536)^ mass is also presented. 

These measurements are performed by fitting the invariant mass spectrums of combinations ® of 
D^TT^, Dfj, Df-K^j, and Df'K'^TT~ particles. Combinations of Dfir'^ and Dfir' are also studied in 
the search for new states. The following two sections of this paper describe those details in common 
to the entire analysis. The study of each combination is then discussed in their own section. The 
paper ends with a summary. 

2 CANDIDATE SELECTION 

This analysis is performed using a 125 fb^^ data sample collected on or near the T(45) resonance 
with the BABMi detector at the PEP-II asymmetric-energy e"'"e~ storage rings. The BABAR detector, 
a general-purpose, solenoidal, magnetic spectrometer, is described in detail elsewhere [6]. Charged 
particles were detected and their momenta measured by a combination of a drift chamber (DCH) 
and silicon vertex tracker (SVT), both operating within a 1.5-T solenoidal magnetic field. A 
ring-imaging Cherenkov detector (DIRC) is used for charged-particle identification. Photons are 
detected and measured with a Csl electromagnetic calorimeter (EMC). 

All of the final states explored in this analysis involve one Df candidate decaying to K~^K~'k~^ . 
A clean sample of candidates is obtained using particle identification by requiring a Cherenkov 
photon yield and angle consistent with the hypothesis. This information is augmented with 
energy loss measurements in the tracking systems. The efficiency of identification is approx- 
imately 85% in the kinematic range used in this analysis with a vr"^ misidentification rate of less 
than 2%. A similar procedure is used to produce a sample of vr^ candidates. 

Each Df candidate is constructed by combining a K~^K~ candidate pair with a tt"*" candidate 
in a geometrical fit to a common vertex. An acceptable K'^K~'k'^ candidate must have a fit 
probability greater than 0.1% and a trajectory consistent with originating from the e^e~ luminous 
region. Backgrounds arc further suppressed by selecting decays to K*^K^ and ^vr^. Additional 
details of this selection procedure can be found elsewhere [1]. Combinations of K'K^-k^ with 
1.954 < m{K^K-'K+) < 1.981 GeV/c^ are taken as D+ candidates. 

For the purposes of calculating the invariant mass of the various particle combinations in this 
paper, the energy of each Df candidate is calculated from the measured momenta and the PDG 
Df mass of 1968.5 MeV/c^ [7]. The uncertainty in this mass (0.6 MeV/c^) is taken as a systematic 
uncertainty. 

Once a Df candidate is obtained, a search is performed for all accompanying 7r°, 7, and tt^ 
particles. All energy clusters in the EMC unassociated with charged tracks and consistent with an 

^Inclusion of charge conjugate states is implied throughout this paper. 
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electromagnetic shower are considered photon candidates. A candidate vr^ is formed by constraining 
a pair of photons each with energy greater than 100 MeV to emanate from the intersection of the 
D'^ trajectory with the beam envelope, performing a one-constraint fit to the 7r° mass, and requiring 
a fit probability greater than 2%. 

A given event may yield several acceptable tt" candidates. For the various final-state samples 
used in this analysis, we retain only those candidates for which neither photon belongs to another 
otherwise acceptable 7r° of momentum greater than 150 MeV/c. In a similar fashion, we discard 
any 7 candidate which belongs to the same list of acceptable 7r° candidates. 

3 SIMULATION AND CALIBRATION 

The efficiency and reconstructed mass distribution of all final-state combinations discussed in this 
document are derived from detailed simulation based on the Geant4 [8] toolkit, reconstructed with 
the same algorithms as the data, and passing the same selection requirements. Individual samples 
of between 120,000 and 240,000 events were generated for each individual signal. 

In addition, for the purposes of optimizing candidate selection requirements, a simulated sample 
corresponding to 30 fb~^ of e'^e" cc events is used to describe combinatorial background. 

Clean samples of r — > pi^ and t ^ ttu were studied to validate the accuracy of the simulation of 
the reconstruction efficiency and energy scale for photon and -k^ candidates. Based on this study, a 
photon efficiency correction of 1.6% and EMC energy scale corrections of between 0.2% and 0.8% is 
applied. The remaining systematic uncertainty in the Monte Carlo prediction of neutral efficiency 
and energy scale is a relative 3% and 1%, respectively. 

The reconstruction and selection efficiency of the various final-state combinations studied in 
this analysis tend to be dependent on the center-of-mass momentum p* of the D*j(2317)"'" and 
Dsj{2AQQ)'^ . To minimize any dependence on the p* distribution assumed in the Monte Carlo, 
each final-state is restricted to the same minimum p* value of 3.2 GeV/c^ (this also serves to 
remove production from B decay). In addition, a p* dependent efficiency function derived from 
Monte Carlo is used to calculate efficiency corrected yields for all decay modes. 

Estimates of the tracking efficiency obtained from r decays and from a comparison of SVT and 
DCH measurements indicate that the Monte Carlo prediction of tt^ efficiency is, on average, 0.8% 
too high, with a systematic uncertainty of 1.3%. 

4 THE D+TT^ SYSTEM 

To form D^tt^ combinations, each Df candidate is combined with one ir^ candidate. The p* of the 
resulting combination is required to be greater than 3.2 GeV/c. The corresponding invariant mass 
distribution is shown in Fig. 1. Signals for D*j(2317)+ and Z)*(2112)+ are clearly visible. 

For the purposes of fitting the D^n^ mass spectrum and extracting information on the -D*j(2317)' 
and -Dsj(2460)^ mesons, the Dfir^ sample is restricted to those vr*^ candidates with momentum 
greater than 400 MeV/c. This momentum requirement is based on an optimization of Q = ej^fB^ 
where e is the efficiency for the hypothetical decay Dsj(2460)"'" Dfir'^, and B is the amount of 
background as predicted by Monte Carlo. The resulting invariant mass distribution is shown in 
Fig. 2. As a result of this stricter tt^ momentum requirement, the L'*(2112)"'" — > D^tt^ signal is 
entirely eliminated. 

Estimates of the yield and mass of the D*j(2317)+ are obtained using a unbinned likelihood 
fit of the Df-K^ mass spectrum. The lineshape is also used to calculate a limit on the yield from 
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Figure 1: The D^-k^ invariant mass distribution for candidates with loose 7r° requirements. 
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Figure 2: The Dfn invariant mass distribution for candidates that satisfy the requirements dis- 
cussed in the text. The solid curve is the result of a unbinned likelihood fit. The dark (light) gray 
region is the contribution from the J (2460)+ Dfir^^ {D*{2112)+ Df-f) reflection. The 
inset is an expanded view near the D*j(2317)+ mass. 
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Figure 3: The kinematics of Ds.j{2A&{))^ D'^tt^^ decay. The shaded area indicates the region of 
Dfir^ and Df'y mass that is kinematically allowed in the decay of an object of mass 2459.2 MeV/c^. 
The lines mark the kinematic space limited to the decay which proceeds through intermediate 
£)*(2112)+ or D*j{2317)+ mesons. 



hypothetical Dsj{2460)~^ Dfir^ decays. The lineshapcs of the reconstructed D*j{2317)~^ and 
L'sj(2460)^ mass distributions are obtained from fits to simulated signal events. In both the 
simulation and in the fit it is assumed that the intrinsic width T of both mesons is small enough 
that this lineshape is not significantly affected. 

The background to the D*j{23n)~^ comes from unrelated Df and tt^ (the combinatorial back- 
ground) and from two types of reflections. One reflection, from D*{2112)'^ Dfj decays com- 
bined with an unassociated 7 to form a fake vr'' candidate, is peaked near the kinematic limit in 
Df-K^ mass of 2154.6 MeV/c^. The second reflection is produced by the D'^ and 7r° mesons from 
Dsj(2460)+ — ^ I?*(2112)"'"7r'' decay. Due to a kinematic coincidence, this reflection has a mean 
mass that is close to the 15*^(2317)^ mass and must be accurately determined in order to correctly 
measure the D*j(2317)+ properties. As shown in Fig. 3, if reconstruction efficiency is uniform, 
the Dsj{2A&0)^ reflection would appear as a nearly square function of total width 41.5 MeV/c^ 
smeared by resolution. Because of the variation in reconstruction efficiency across the Dsj(2460)"'" 
phase space, the reflection shape is asymmetric. The prediction from Dsj{24QQi)^ Monte Carlo is 
used to determine the precise shape. 

The likelihood fit to the Dfi:^ mass spectrum is shown in Fig. 2. Whereas the shape and 
magnitude of the Dsj(2460)"^ reflection are fixed to Monte Carlo predictions, the parameters of 
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the D*(2112)"'" reflection and the shape of the combinatorial background are allowed to vary in 
the fit. The mass and yield of the D*j(2317)+ are determined by shifting the associated lineshape 
up or down in mass and adjusting the overall amplitude to best match the data. A limit on the 
yield of the hypothetical Dsj{2AQ0)~^ Dfir^ decay is determined by adjusting the amplitude of 
the associated lineshape positive or negative to best match the data. The resulting likelihood fit 
successfully describes the data. 

The fit determines (statistical errors only) a D*j{2317)+ mass of [2318.9 ±0.3] MeV/c^ and raw 
L>*j(2317)+ and I?,,,/ (2460)+ yields of 1275 ± 45 and 3 ± 26 mesons. 

A systematic uncertainty specific to the fit to the -D+tt" system is the size of the Dsj(2460)"'" 
reflection. If the likelihood fit is allowed to adjust the size of this reflection to best match the data, 
the result is a D*j(2317)+ yield 0.5% smaller and a shift in D*j(2317)+ mass that is less than 
0.1 MeV/c^. In addition, various different models of the L'*(2112)"'" reflection can be used in the 
fit, in which case variations in the D*j(2317)"'" yield of up to 0.4% are observed. 

5 THE D+-f SYSTEM 

To form combinations, each D+ candidate is combined with 7 candidates with energy greater 
than 500 MeV/ c. This energy requirement is based on an optimization of Q' = S/ \/5"+~B, where 
S is the expected signal size [3] and B is the amount of background as predicted by Monte Carlo 
for the decay Dsj{2460)'^ — Dfj. In addition, any Df^ combination with p* less than 3.2 GeV/c 
is discarded. The resulting invariant mass distribution is shown in Fig. 4. 

For reasons of simplicity, the likelihood fit is restricted to a -D+7 mass range between 2.15 
and 2.85 MeV/c^. As shown in Fig. 4, there are two structures clearly visible in this spectrum 
within this range on top of a gradually falling background distribution. The higher mass structure 
corresponds to the Dsj(2460)"'" meson. The likelihood fit uses a signal lineshape determined from 
signal Monte Carlo. To determine the Z)sj(2460)+ mass, the mean value of the signal is allowed 
to shift. To determine the Dsj(2460)+ yield, the fit determines the amplitude that bests describes 
the data. 

The lower mass structure is a combination of two reflections. The largest reflection is composed 
of the meson from the decay Z)*j(2317)+ Df-K^ combined with one of the photons from the 
tt". The second is produced in a similar fashion from Z)sj(2460)+ Dfn^^ decay. The shapes 
of these two reflections are influenced by the candidate selection requirements and, in particular, 
the requirement that the 7 particle not belong to a vr*^ candidate. This interdependency is complex 
enough that the Monte Carlo simulation fails to adequately reproduce the precise shape of the 
reflection. For this reason a parameterization of D*j(2317)+ —>■ Dfir^ reflection shape is included 
as part of the fit. The Dsj{2A60)'^ — iD+7r°7 reflection shape is fixed to the prediction from Monte 
Carlo. 

The kinematics of L'*j(2317)"'" Dfir^ decay require that the Dfj reflection end at Dfj 
mass of 2.304 MeV/c^. The lineshape assumed by the likelihood fit respects this limit but includes 
detector resolution effects that will tend to blur this boundary. Just adjacent to the boundary is 
where a hypothetical D*j(2317)"'" Dfj signal would reside. The lineshape of the (2317)+ 
signal is extracted from signal Monte Carlo. Because of uncertainties in the precise shape of the 
L>*j(2317)+ — > Dfn^ and L>sj(2460)+ Dfn^^ reflections, there is considerable uncertainty in 
the potential amount of Z)*j(2317)+ — -D+7 decays. 

The likelihood fit represents the background by a constant plus the tail of a Gaussian distribu- 
tion. The result of the fit is shown in Fig. 4. A Dsj{2460)+ mass of [2457.2 ± 1.6] MeV/c^ and 
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Figure 4: (top) The Df^ invariant mass distribution for candidates that satisfy the requirements 
discussed in the text. The sohd curve is the result of a unbinned Hkclihood fit. (bottom) The 
mass distribution after subtracting the contribution from combinatorial background as estimated 
by the fit. Various contributions to the signal and reflection portions of the fit are overlaid. The 
-D*j (23 17)+ signal (dot-dash line) appears as a negative contribution due to a negative fluctuation 
in the fit. 
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yield of 509 ± 46 mesons is obtained (statistical errors only). The fit, which allows the signal yield 
to fluctuate to negative values, obtains — 107± 84 D*j(2317)+ decays (statistical errors only). 

The largest uncertainty in the likelihood fit is the shape of the reflection and background. Any 
tails in the signal or reflection mass distributions will tend to fill in the region of mass around 
2.35 MeV/c^, suppressing the size of the background and inflating signal yields. In addition, as 
mentioned above, the D*j(2317)+ yield is sensitive to the shape of the reflections because the two 
peaks overlap. In particular, if the Monte Carlo predictions for the reflect shape and size is used 
in the fit, the D*j{2317)~^ yield is increased to 209 mesons. 

Another systematic check is to use alternate functions to describe the combinatorial background, 
in which case it is possible to change the Dsj{2460)^ yield by at most 3% without adversely affecting 
the fit quality. 

6 THE i:>+7r07 SYSTEM 

A likelihood fit to the Dfir^^ invariant mass distribution is used to determine the mass and yield 
of the £)5j(2460)+ meson. Because this is a three body decay, there are additional issues that need 

to be addressed. In particular, it is reasonable to assume that the Dsj{2A60)+ decays to D+tt^j 
through either of two possible intermediate resonances: D*{2112)^-k^ or L'*j(2317)^7. As shown 
in Fig. 3, these two decay modes overlap in all three mass projections. There is no a priori reason 
to favor one decay path over another. A second, two-dimensional fit is used to distinguish them. 

To select the data sample for studying the -Dsj(2460)"'" meson, the quantity Q' = S/^/S + B is 
optimized simultaneously for both the minimum vr'' momentum and 7 energy. The expected signal 
size S is calculated based on results from previous measurements [5]. Based on this procedure a 
minimum tt^ momentum of 400 MeV/c and minimum 7 energy of 135 MeV is chosen. 

Included in Fig. 5 is the Z?+7r°7 mass spectrum for all selected candidates below a mass of 
2.75 GeV/c^. A peak near the Dsj(2A60)~^ mass is apparent. The signal shape of the _Dsj(2460)^ 
is obtained by fitting a signal Monte Carlo sample. This shape is combined with a second order 
polynomial to fit the mass spectrum. The result is a Dsj{2460)'^ mass of [2457.8 ±2.8] MeV/c^ and 
yield of 246 it 44 mesons (statistical errors only). The Z?sj(2460)"'" peak, however, is underlayed by 
a substantial background which has structure that is not well represented by a simple polynomial. 
To obtain a more accurate measurement of the -D^j (2460)+ meson requires a more refined selection 
to isolate the signal. 

Since both the Z)*(2112)+7r° and I?*j(2317)+7 decay modes are restricted to a Dfj mass near 
the D*(2112)+ mass, the Dsj{24:60)~^ signal can be isolated by requiring the invariant mass 
to reside within 2 MeV/c^ of the D*(2112)+ mass (2112.4 MeV/c^ [7]). The result is included in 
Fig. 5. Although the background is now substantially reduced, it is clear that some type of peaking 
background has been introduced by this requirement since the size of the peak in the restricted 
data sample is noticeably larger than that of the inclusive sample. 

The source of the peaking background can be investigated further by considering two Z)*(2112)"'" 
sidebands, each of total width 40 MeV/c^ in Df^ mass and with centers separated by ±60 MeV/c^ 
from the D*(2112)"'" mass. The Dfir^^ mass distribution in these sidebands is included in Fig. 6. 
Two reflections are observed. Both can be identified in the Dfir^ invariant mass distribution of 
Fig. 2. The peak that appears near a Dfn^^ mass of 2.25 MeV/c^ in the lower -Dj"7 sideband 
originates from £)*(2112)+ — >■ Dfj decays combined with two unassociated 7 particles, one of 
which has been combined with the 7 from the D*{2112)~^ decay to form a fake tt^ candidate. A 
calculation of kinematics based only on the Df, 7r°, and (2460)"*" masses can demonstrate that 
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Figure 5: The inclusive and semi-inclusive Dfir^'j invariant mass distributions. Shown in solid 
points is the total data sample. Overlaid is the simple line fit discussed in the text. Shown in open 
round points is the subset which has a Df'j mass within ±2 MeV/c^ of the Z)*(2112)+ mass. 
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the D*{2112)~^ reflection will appear with a predictable amount of smearing at a specific Dfir^'j 
mass in the Dfj signal and sideband selections. 

A more important reflection is from D*j(2317)''" ^^7r° decays combined with unassociated 
7 particles. The position and lincshape of this reflection has been determined using D*j{2317)~^ 
Monte Carlo. Although the Monte Carlo prediction of the size of the reflection is consistent with 
the data within statistics, its size has been adjusted uniformly upwards by 23% (corresponding 
to approximately one standard deviation) to best match the amplitudes observed in the two Dfj 
sidebands. The accuracy of the Monte Carlo prediction of the D*j(2317)'^ shape has already been 
confirmed in the Dfn^ fit shown in Fig. 2 

Independent likelihood fits are applied to each of the sideband distributions and to the signal 
distribution. All three fits can successfully model their respective data samples, as shown in Fig. 6. 
The fit to the D^^ signal region includes those two reflections and the Dsj{24:60)~^ signal. The 
Dsj{2460)'^ mass is included in the fit by allowing the lineshape to up or down in mass. The fit 
obtains a L>sj(2460)+ mass of [2459.1± 1.3] MeV/c^ and yield of 292±29 mesons (statistical errors 
only). 

A systematic uncertainty specific to this fit is the size of the D*j(2317)"'" reflection. Using the 
prediction of yield from the Monte Carlo unchanged increases the L'sj(2460)"'" yield to 308 mesons 
and shifts the mass upwards by 0.3 MeV/c^. 

The fit of the Dfir^j mass does little to differentiate between the two possible Dsj(2460)"^ 
decay modes L'*(2112)"'"7r'' and D*j{2317)~^^. It is necessary for this purpose to investigate the 
L>+7r° and Df 7 mass distributions more closely. As shown in Fig. 3, the D* {2112)~^ tt^ decay can be 
expected to appear as a sharp peak in the Df^ spectrum centered on the L>*(2112)+ mass, and with 
a distribution spread out by approximately 41.3 MeV/c^ in Dfir^ mass centered at 2313.4 MeV/c^. 
In contrast, the decay -Dsj(2460)"'" — > I?*j(2317)"'"7 appears as a sharp peak in the Dfn^ spectrum 
and as a distribution spead out by approximately 40.9 MeV/c^ in the Dfj spectrum centered at 
2110.5 MeV/c^. 

The Dfir'^ and Df^ mass distributions of the signal cannot be explored without correctly sub- 
tracting backgrounds from D*(2112)+ L>+7 and L>*j(2317)"'" Dfn^ decays. This subtraction 
is performed by a two-dimensional unbinned likelihood applied to the Dfir^ and Df^ mass distri- 
butions of the data. The likelihood fit is restricted to the data sample contained in the boundaries 
illustrated in Fig. 7. This fit includes five sources of Dfir'^j candidates: 

• Combinatorial background represented by a two-dimensional quadratic function. 

• Background from L'*(2112)^ Dfj decays combined with unassociated vr*^ candidates rep- 
resented by a L'*(2112)+ lineshape in the Df^ mass and as a linear function in Dfir^ mass. 

• Background from D*j(2317)"'' Dfii^ decays combined with unassociated 7 candidates 
represented by a -D*j(2317)''~ lineshape in the Dfir^ mass and as a linear function in Dfj 
mass. 

• A signal from Dsj{2A60y L»*(2112)+7r° represented by a D*(2112)+ lineshape in the Dfj 
mass and a 41.5 MeV/c^ wide reflection in the Dfir'^ mass with a specific shape and mass 
predicted by signal Monte Carlo. 

• A signal from Dsj(2460)+ D*j(2317)+7 represented by a L>*j(2317)+ lineshape in the 

DfiT^ mass and a 40.9 MeV/c^ wide reflection in the Dfj mass with a shape determined 
from kinematics smeared by a resolution extracted from signal Monte Carlo. 



17 



O 50 





EiEiR 

Preliminary 1 25 fb 

' JmJ- 


i 











2.1 2.2 2.3 2.4 2.5 2.6 2.7 

D* 7T° 7 Invariant Mass (GeV/ c^) 

Figure 6: The Dfir^'y mass spectrum of the sample of events that fall in a Df'y signal region 
(center plot), the Dfj high mass sideband (top plot), and the low mass sideband (bottom plot). 
The -057(2460)"*" signal appears only in the signal region. The results of separate likelihood fits are 
superimposed on all three distributions. Reflections from (dark gray) Z)*j(2317)^ Dfir^ and 
(light gray) Z)*(2112)+ decays appear in all three distributions at positions that are fixed 

from kinematics. The relative contribution of the D*j(2317)"'" reflection in each sample has been 
determined using signal Monte Carlo. 
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Figure 7: A scatter plot of the Dfir^ versus Dfj invariant masses from the sample of Dfn^^ 
candidates. Two bands from D*j(2317)+ D+tt^ and L>*(2112)+ Dfj decays are clearly 
visible as background. The D5j(2460)^" signal consists of an excess of candidates near the area 
these two bands cross. The grid lines indicate the range of events used in a likelihood fit. 
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The results of this Ukehhood fit is shown in Fig. 8, divided into the regions dehneated by the 
grid shown in Fig. 7. The fit produces an adequate model of the data in all regions. The result 
(statistical errors only) is a r>;(2112)+7r° yield of 266 ± 38 mesons and a L»*j(2317)+7 yield of 
— 11 lb 37 mesons, the former number being somewhat smaller than the yield determined by the 
first fit, though consistent within systematic uncertainties. Based on these results, it appears that 
the decay Z'sj(2460)+ Dfir^j can be successfully described as proceeding entirely through the 
channel i:»*(2112)+7r0. 

The likelihood fit of Fig. 8 is sensitive to various assumptions. For example, the D*j(2317)~''7 
yield can be raised to 2 candidates while preserving the quality of the fit if the mass resolutions 
are degraded slightly. In addition, possible inaccuracies in the inclusive D*(2112)+ p* distribution 
predicted by Monte Carlo can produce changes in the £)*(2112)+ lineshape that can raise the yield 
of D*(2112)+7r° decays to 280. 

7 THE D+7T+7r- SYSTEM 

To form £)+7r+7r~ candidates, each Df is combined with 7r+ and 7r~ candidates with momentum 
above 250 MeV/c. This momentum requirement is obtained from the optimization of Q for the 
decay Dsj{2460)~^ D+tt+tt". The resulting invariant mass distribution has two distinct, narrow 
peaks, as shown in Fig. 9. These two peaks correspond to the decays of the Dsj{2460)~^ and 
£)si(2536)+ mesons, as previously noted by the BELLE collaboration [3]. 

There is a type of reflection that could be introduced into the inclusive distribution of Z^^vr^vr" 
mass that would not appear in the other combinations discussed up to this point. An example 
of such a reflection is the hypothetical decay of a heavy particle into or K* plus two charged 
particles, either of which might he a K meson or proton. The corresponding (j) or K* from this 
decay can be combined with an unassociated tt~^ candidate to form a flctitious D'^ candidate. The 
result could produce a peak in the Dfir^ir^ mass distribution when combined with the other decay 
products. As a sanity check, the Df decay products are combined with the tt^ candidates using 
various different particle species hypotheses (tt^, p, K^) to check for underlying resonances not 
associated with Df decay. No such resonances are uncovered. 

The likelihood fit to the DfTT^TT~ mass distribution consists of three signal distributions 
(Z)*j(2317)+, Z)sj(2460)+, or Dsi(2536)+ decay) plus a third-order polynomial to describe the 
combinatorial background. The shapes of the three signals are derived from signal Monte Carlo 
samples. The Dsj(2460)~^ and Dsi{2536)~^ signal shapes are allowed to shift upwards and down- 
wards in mass in order to best describe the data. The result of this fit, included in Fig. 9, is 
(statistical errors only) a D*j(2317)"'" yield of 0.6 ± 1.8 decays; a Dsj(2460)^ mass and yield of 
[2460.1±0.3] MeV/c2 and 67±11 decays; and a Dsi(2536)+ mass and yield of [2534.3±0.4] MeV/c^ 
and 124 ± 18 decays. 

A systematic uncertainty specific to this likelihood fit is the assumption of the background 
shape. Substituting a fourth-order polynomial for the background changes the Dgj (2A60)'^ and 
D*(2112)+ yields by a relative 5.3% and 9.1%. 

8 THE D+TT"^ SYSTEM 

There has been some conjecture [9, 10, 11] that the D*j{2317)~^ may be a four-quark hybrid state. 
It might be expected, if this was true, that neutral and doubly-charged partners should exist with 
a similar mass. The Dfir^ system can be used to test this possibility. 
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Figure 8: The results of a likelihood fit to the Dfir^ and Dj"7 masses of the Dfir^j candidates 
which appear inside the grid shown in Fig. 7. Each plot corresponds to a column or row of that grid 
as indicated by the numbers in parentheses. The points correspond to the data. The histograms 
are a Monte Carlo representation of the fit results. Shown in gray (black) are the contributions 
from Dsj{2A60)+ i:>|(2112)+7r0 (£),j(2460)+ ^ L»*j(2317)+7) decay predicted by the fit. The 
light gray area is the contribution from combinatorial background. 
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Figure 9: The Df'K~^'K~ invariant mass distribution for candidates that satisfy the requirements 
discussed in the text. The soUd curve is the result of a unbinned hkehhood fit. 



To form Dfir^ combinations, each candidate is combined with vr^ candidates with momen- 
tum greater than 300 MeV/c. This momentum requirement is obtained from the optimization of 
Q for the hypothetical decay D*j{22>n)^ — > Dfir'. The resulting mass distributions are shown in 
Fig. 10. No resonant structure is observed. 

Likelihood fits have been performed assuming a hypothetical D*j(2317)^ or D*j{23n)~^~^ state 
at precisely the L'*j(2317)+ mass. The results of these fits are included in Fig. 10. The yields 



obtained for the fit (statistical errors only) for the D*j(2317)'' and D*j(2317) 
—39 ± 16 mesons, respectively. 



are -28 ± 25 and 



9 SYSTEMATIC STUDIES 

There are various systematic uncertainties in common to many of the likelihood fits presented 
here. As discussed earlier, there is an uncertainty of ±0.6 MeV/c^ in the given Df mass. A 
1% uncertainty in the energy calibration of the EMC results in uncertainties of between 0.6 and 
0.8 MeV/c^ for those combinations which include a vr^ or 7 candidate. The largest uncertainty in 
the vr^ momentum arises from the accuracy of the calculations applied by the track reconstruction 
software to account for the energy loss in the material of the detector. A conservative estimate of 
the effect on the I?+7r+7r~ mass of all tracking uncertainties is ±1 MeV/c^. 

The mass lineshapes used in the likelihood fits of this analysis, derived from Monte Carlo 
simulation, provide models of the data of satisfactory quality. Variations in the resolution of the 
detector have been applied to widen these lineshapes within reasonable bounds. The resulting 
change in yield or mass is assigned as a systematic uncertainty. 

As mentioned earlier, the reconstruction efficiency for all particle combinations discussed in this 
paper vary as a function of p*. Yields are calculated by performing separate likelihood fits with 
candidates weighted by the inverse of this efficiency. To check this approach, likelihood fits are also 
performed on the candidate sample divided into bins of p*. The yields from the binned data are 
calculated using the average reconstruction efficiency in that bin. Any difference in yield produced 
by these two methods is assigned as a systematic uncertainty. 
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Figure 10: The DfTr~ (top) and Dfn'^ (bottom) invariant mass distributions for candidates that 
satisfy the requirements discussed in the text. The solid curve is the result of a unbinned hkelihood 
fit. 
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Any inaccuracies in the Monte Carlo predictions of the reconstruction efficiency will affect 
the calculation of branching ratios. As already discussed, the uncertainty in 7r° or 7 efficiency is 
assumed to be 3% and the uncertainty in tracking efficiency is assumed to be 1.3%. 

10 COMBINED RESULTS 

A preliminary measurement of the D*j(2317)"'" mass has been obtained from the Dfi:^ system: 

m(D;j(2317)+) = 2318.9 ± 0.3 (stat.) ± 0.9 (syst.) MeV/c^ . (1) 

Preliminary measurements of the Dsj(2460)+ has been obtained from Df^, Dfn'^^, and 
-D^TT'^Tr" decays. The results, in units of MeV/c^, are: 



m(Dsj(2460)+ ^ D+j) = 2457.2 ± 1.6 (stat.) ± 1.3 (syst.) (2) 
m{Dsj{2A60)+ Dfir^j) = 2459.1 ± 1.3 (stat.) ± 1.2 (syst.) (3) 
m(Dsj(2460)+ ^ D+Tr+Tr") = 2460.1 ± 0.3 (stat.) ± 1.2 (syst.) . (4) 



The average of these results, after accounting for correlated systematic uncertainties, is: 

m{Dsj{2m)+) = 2459.4 ± 0.3 (stat.) ± 1.0 (syst.) MeV/c^ . (5) 

The fit to the -D+tt+tt"" mass spectrum also provides a preliminary estimate of the Z)si(2536)'*" 
mass: 

m(Lisi(2536)+) = 2534.3 ± 0.4 (stat.) ± 1.2 (syst.) MeV/c^ . (6) 

The following preliminary branching ratios have been estimated by comparing the efficiency 
corrected yields from the corresponding decay modes. 

Br(aX2460)+ Dfy) _ ^^^^^ ^ ^^^^^^ ^ ^^^^^ ^^^^^ p, 



Br (L>, J (2460)+ Dtn'^-f) 

MDsjiMeor^Dt.Y) ^ ^ ^ ^^^^ 

Br{Dsj{2460)+ D+n^j) ^ ^ ^ ^ ^ 

A 95% CL upper limit has been calculated for the decay modes in which no significant signal is 
observed. This limit is based on the deviation from the measured yields of 1.96 standard deviations, 
derived from the quadrature sum of statistical and systematic uncertainties. 



l3r{D*j{2317)+ ^ D+ttO) 

Br{Dsj{2A60)+ D+t:^) 
Br(L>,j(2460)+ ^ DtnO^) 

Br{Dsj{2m)+ ^ £>;j(2317)+7) 
i3r(L',j(2460)+ ^ DfTr^-f) 



< 0.11 (10) 
< 0.23 (11) 



Br{Dtj{2317)+ ^ D+tt+tt-) , 
BriD*j{2317)+ ^ D+tt^) 



(12) 
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11 CONCLUSION 



Based on 125 fb"^ of e+e" cc data, the i:>;j(2317)+ and Dsj(2460)+ masses are measured to be 
[2318.9±0.3 (stat.)±0.9 (syst.)] MeV/c^ and [2459.4±0.3 (stat.)±1.0 (syst.)] MeV/c^, respectively. 
Both values are consistent with previous measurements [2, 3, 4] with perhaps the exception of the 
£)*j(2317)+ mass which is 1.6 MeV/c^ higher than our previous result [5]. 

Significant signals are observed for L>*j(2317)+ ^ D+tt^, L»^j(2460)+ ^ Df-f, L»sj(2460)+ ^ 
Z)+7r°7, and Dsj{2AQ0)+ D+Tr+Tr". The data is consistent with the Aj(2460)+ D+n^j 
decay proceeding entirely through the channel D*(2112)"'"7r'^. The Dfir^TT^ result confirms the 
decay first observed by Belle [3], although at a rate that is somewhat lower. 

No significant signal is observed for cither the Df^ or DfTr'^Tr~ decays of the D*j(2317)'''. A 
search for the Dfir^ decay of the Dsj(2460)+ also produces no significant signal. Searches for a 
neutral or doubly-charged partner for the D*j{2317)^ decaying to Dfn^ or Dfir' also produce 
no significant signal. 
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